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ABSTRACT 

The  stability  of  two-dimensional  infinitesimal  disturbances  of  the 
inviscid  Karman  vortex  street  of  finite  area  vortices  is  reexamined. 

Numerical  results  are  obtained  for  the  growth  rate  and  oscillation  frequencies 
of  disturbances  of  arbitrary  subharmonic  wavenumber  and  the  stability 
boundaries  are  calculated.  The  stabilization  of  the  pairing  instability  by 
finite  area  demonstrated  by  Saffman  and  Schatzman  (1982)  is  confirmed  and  also 
Kida's  (1982)  result  that  this  is  not  the  most  unstable  disturbance  when  the 
area  is  finite.  Contrary,  however,  to  Kida's  quantitative  predictions,  it  is 
now  found  that  finite  area  does  not  stabilize  the  street  to  infinitesimal  two¬ 


dimensional  disturbances  of  arbitrary  wavelength  and  that  it  is  always 
unstable  except  for  one  isolated  value  of  the  aspect  ratio  which  depends  upon 
the  size  of  the  vortices.  ■***■• 
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THE  STABILITY  OF  INVISCID  VORTEX  STREETS  OF  FINITE  CORED  VORTICES 
D.  I.  Meiron*,  P.  G.  Saffman**,  and  J.  C.  Schatzman*** 

1.  Introduction 

We  consider  the  stability  to  two-dimensional  Infinitesimal  disturbances 
of  the  steady  lnvlscid  incompressible  flow  produced  by  Infinite  rows  of  finite 
cored  vortices.  The  undisturbed  flow  In  a  frame  of  reference  moving  with  the 
vorticity  Is  described  by  a  stream  function  ¥(x,y)  with  the  property 
that  ¥(x  +  i  ,y)  5  ¥(x,y),  and  ¥(x,“)  "  U^y,  ¥(x,-«)  ■  U^y.  *r*  two 

cases:  (1)  wake  type  flows  In  which  ^1  "  ^2  "  8*y»  an^  (11)  mixing  layer 

type  flows  with  U2  -  -Oj  ■  0.  In  the  first  case,  there  Is  no  net  vorticity; 
the  canonical  flow  Is  that  of  the  Karman  vortex  steeet  of  two  Infinite 
straight  parallel  staggered  rows  of  point  vortices  of  equal  and  opposite 
circulation  T ,  with  separation  1  parallel  to  the  rows,  distance  h  between 
the  rows,  and  each  vortex  opposite  the  midpoint  of  vortices  In  the  other 
row.  In  the  second  case,  there  Is  net  vorticity  and  the  canonical  example  is 
a  single  straight  infinite  row  of  point  vortices  of  equal  circulation  T 
separated  by  distance  1  . 

Infinitesimal  disturbances  to  these  flows  are  described  by  a  perturbation 
cT'(E»y>c)  to  the  stream  function,  where  e  is  Infinitesimal  and  by  Floquet  or 
Bloch  wave  theory  the  perturbation  is  a  sum  of  modes  of  the  form 

j  ,  (,).*! ««  .  (1.i, 


If  follows  from  the  Euler  equations  that  the  perturbation  satisfies  the  linear 


partial  differential  equation 


0  . 


(1.2) 


Substitution  of  (1.1)  into  (1.2)  and  requiring  that  T '  be  bounded  gives  an 
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eigenvalue  problem  for  the  unknown  growth  rate  o  as  a  function  of  the 
undisturbed  flow  and  p.  Notice  that  p  la  an  arbitrary  real  number,  which 
clearly  without  loss  of  generality  may  be  taken  In  the 

range  0  <  p  <  1  or  p  <  ^  »  ®l°ce  adding  an  integer  to  p  Is  simply 

equivalent  to  relabeling  the  eigenvector  If  p  ■  0,  the  disturbance  has 

the  setae  spatial  period  1  as  the  undisturbed  flow,  and  will  be  called  a 
superharmonic  disturbance.  If  p  *  0,  the  disturbance  has  in  general  a 
component  of  wavelength  i/p  or  i/(l~p)»  which  Is  larger  than  1,  and  will 
be  called  subharmonic.  The  case  p  ■  1/2,  In  which  the  disturbance 
wavelength  Is  double  that  of  the  undisturbed  flow,  has  attracted  particular 
attention  In  flows  of  mixing  layer  type  where  it  Is  called  the  pairing 
Inatablllty . 

If  a  la  pure  Imaginary  for  all  values  of  p,  the  disturbance  does  not 
grow  with  time  and  the  flow  la  said  to  be  stable.  (Strictly  speaking,  It  Is 
only  stable  to  Infinitesimal  disturbances,  but  the  behavior  of  finite 
amplitude  disturbances  Is  an  open  and  difficult  question  and  will  not  be 
considered  here.)  If  a  has  a  positive  real  part  for  some  value  of  p,  then 
the  flow  Is  unstable. 

The  case  of  point  vortices  has  been  Investigated  in  detail  by  many 
workers;  Lord  Kelvin  and  J.  J.  Thomson  studied  the  single  row  In  the  form  of  a 
circle  (see  Havelock  (1931)  for  a  complete  treatment).  Karman  worked  out  the 
stability  of  the  double  row;  Lamb  (1932)  provides  an  account  of  this  work  and 
describes  the  stability  of  the  single  row  which  Is  always  unstable.  The 
staggered  double  row  Is  likewise  unstable  except  when 
h/i  ■  *c  ■  cosh  *(/2)/ir  ■  0.280550,  at  which  value  Re(o)  ■  0. 

Recently,  due  partly  to  the  availability  of  large  scale  computing 
resources,  the  case  in  which  the  vortices  are  of  finite  size  has  attracted 


attention.  Saffaan  and  Santo  (1981)  examined  the  single  row  of  uniform  finite 
else  vortices,  where  the  vortlcity  Is  constant  Inside  each  vortex  and  aero 
outalde.  The  stream  function  Is  continuous  and  has  continuous  first 
derivatives,  but  the  second  derivatives  have  simple  jumps  across  the 
boundaries  of  the  vortices.  They  calculated  steady  shapes  numerically,  and 
discussed  the  stability  to  two-dimensional  disturbances  by  global  energy 
methods,  which  have  the  advantage  that  they  are  net  limited  to  infinitesimal 
disturbances.  They  concluded  that  superharmonic  disturbances  are  stable  when 
the  vortices  are  not  too  large,  but  there  exists  a  critical  sise  at  which 
superharmonic  disturbances  become  unstable.  This  provides  another  mechanism 
for  the  evolution  of  a  single  row  which  goes  under  the  name  of  'tearing' .  The 
analysis  showed  that  the  pairing  Instability  would  not  be  eliminated  by 
effects  of  finite  else,  but  the  qualitative  global  analysis  could  not 
determine  whether  the  growth  rate  Increased  or  decreased.  (Closed  form 
quantitative  eolutlons  for  the  case  of  hollow  vortices  suggested  that  the 
effect  of  slxe  on  the  pairing  Instability  la  small,  see  also  Baker,  Saffman 
and  Sheffield  (1976)).  Numerical  solutions  by  Plerrehumbert  and  tfldnall 
(1981)  for  the  so-called  Stuart  vortices  show  similar  behavior  (this  paper 
also  contains  results  for  three-dimensional  disturbances  which  are  outside  the 


scope  of  the  present  work) . 

The  effect  of  finite  core  size  on  the  stability  of  the  Karma n  vortex 
street  was  considered  by  Doom  (1955),  but  his  undisturbed  flow  was  neither  an 
exact  solution  of  the  Euler  equations  nor  of  the  Navler-Stokes  equations.  As 
will  be  emphasized  below,  the  question  of  the  stabilization  by  finite  size  is 
rather  delicate,  and  Doan's  conclusion  that  finite  size  does  not  stabilize  Is 
of  uncertain  significance.  The  first  consistent  calculation  was  carried  out 
by  Saffaan  and  Schatzman  (1982) ,  who  calculated  the  values  of  a  for 


Dist 


de 

Avail  and/or 
Special 


superharmonic  disturbances  to  a  double  row  of  uniform  vortices  and  the 
subharmonic  pairing  type  instabilities  for  p  •  1/2.  They  formulated  the 
problem  exactly  in  terms  of  a  non-linear  integro-dif ferentlal  equation  based 
on  the  'water-bag'  method  (see  Deem  and  Zabusky  (1978)),  whose  solutions  give 
the  steady  shapes  and  whose  Frechet  derivative  gives  the  linear  eigenvalue 
problem  for  the  o.  These  equations  were  solved  numerically,  thereby  giving 
approximate  solutions  to  an  exact  formulation.  It  was  found  that 
superharmonic  disturbances  were  stable,  as  can  also  be  predicted  by  a  global 
argument,  and  that  the  pairing  instability  could  be  stabilised  by  finite  size 
for  a  range  of  values  of  the  aspect  ratio  h/1  around  ic^.  Saffman  and 
Schatzman  argued  on  the  basis  of  a  symmetry  argument  that  the  pairing 
disturbances  were  the  most  unstable,  at  least  for  small  vortices,  and  it  was 
therefore  not  necessary  to  consider  p  *  ^  with  regard  to  the  stability  to 
general  disturbances.  Unfortunately,  as  will  be  discussed  below,  this 
conclusion  is  false. 

The  problem  was  considered  independently  by  Klda  (1982)  in  a  rather 

different  way.  He  developed  a  perturbation  expansion  in  the  size  of  the 

vortices  to  obtain  approximate  equations  for  the  motion  of  the  vortex 

centroids.  In  order  to  do  this,  it  is  necessary  to  calculate  the  deformation 

of  the  vortex  cores,  which  respond  quasl-steadily  to  the  motion  of  the 

2 

centroids  on  time  scale  t  /T ,  and  can  oscillate  freely  on  the  much  shorter 
time  scale  1/ut  *  A/T ,  where  A  is  the  core  area  and  u>  is  the  magnitude  of  the 
vortlclty  in  the  core.  A  'coarse  graining'  approximation  is  made  which 
neglects  the  high  frequencies  and  produces  approximate  equations  for  the 
vortex  centroids  which  contain  the  finite  size  correction  terms  of 
order  A^/i*.  The  approximate  equations  were  solved  exactly  for  disturbances 
of  arbitrary  subharmonic  wavenumber  p.  Kida  found  for  the  case  p  -  1/2 


-4- 


that  his  results  sgrssd  closely  with  those  of  Saffaan  end  Schatzaan,  despite 

the  dele  by  the  letter  authors  based  on  an  analysis  of  their  results  that  a 

consistent  calculation  of  the  finite  area  stabilisation 

around  h/i  ■  k  required  retention  in  the  dynaaical  equations  of  tens  of 
c 

order  A*/i®.  Hors  iaportant,  however,  was  Kida's  prediction  that  the  most 

unstable  disturbance  does  not  occur  for  p  ■  1/2  but  for  a  value  of  p 

2  4 

differing  by  order  A  /l  ,  in  contradiction  to  the  Saffaan  and  Schatzaan 
syaaetry  arguaent.  He  did  find,  however,  that  the  street  was  still  stabilized 
by  finite  area. 

This  conflict  led  us  to  reexaaine  the  problea  and  tackle  the  case  of 
general  p.  He  will  present  arguaents  and  results  to  show  that  Kida's  claim 
that  the  aost  unstable  disturbances  of  the  Karaan  vortex  street  is 
for  p  *  tyjis  Indeed  correct  and  that  the  Saffaan  and  Schatzaan  syaaetry 
arguaent  was  flawed,  at  least  for  this  case.  (It  is  correct  for  the  single 
row).  On  the  other  hand,  we  shall  also  discuss  the  accuracy  of  the 
perturbation  expansions  and  the  question  whether  neglecting  terms  of 
order  A  /£  in  the  dynamics  is  consistent.  Our  conclusions  here  are  that  it 
is  valid  for  determining  the  effects  of  finite  area  on  the  oscillation 
frequencies  or  growth  rates  of  stable  or  unstable  nodes,  but  not  for 
determining  the  change  of  stability  properties  at  the  critical  spacing 
ratio  kc>  In  fact,  our  calculations  now  predict,  somewhat  surprisingly,  that 
the  street  is  always  unstable,  at  least  for  small  area,  except  for  one 
particular  value  of  h/i  which  depends  upon  the  area.  In  other  words,  the 
qualitative  stability  properties  of  the  point  vortex  configuration  remains 
true  for  finite  area  except  that  the  aost  unstable  disturbance  for  h/i  not 
equal  to  the  special  value  is  not  the  pairing  instability. 

Our  numerical  method  is  sufficiently  general  that  many  different  flows 


can  be  examined  with  little  extra  work.  We  can  therefore  analyse  other  cases 
of  Interest  such  as  the  single  row,  the  symmetrical  double  row,  the  circular 
row  In  which  the  vortices  are  at  the  vertices  of  a  regular  polygon  (this 
provides  a  model  of  the  effect  of  curvature  on  a  shear  layer),  as  well  as  the 
Karman  vortex  street.  He  can  also  handle  the  case  in  which  the  vortices  in 
the  street  have  different  sizes  in  the  two  rows.  It  is,  however,  shown  in 
Appendix  A  that  the  circulations  of  the  vortices  mist  be  equal  and  opposite, 
so  that  steady  motion  of  a  composite  wake  -  mixing  layer  type  flow  in 
which  U2  *  U.  or  *  -0^  is  not  possible,  as  is  also  suggested  by  the 
numerical  calculations  of  Brinlch,  Boldman  and  Goldstein  (1976).  He  shall, 
however,  limit  the  present  paper  to  the  staggered  double  row,  and  present 
results  for  the  other  cases  in  subsequent  papers. 

2.  Consequences  of  symmetry 

He  now  examine  some  of  the  general  properties  that  can  be  established 
Independently  of  the  details  of  the  flow  and  follow  from  the  structure  of  the 
Euler  equations  and  symmetry  of  the  flow. 

Because  equation  (1.1)  is  real,  it  follows  that  if  ¥ '  is  a  mode,  so  is 
its  complex  conjugate  ¥’*.  Thus  if  a  is  an  eigenvalue,  for  a  value  p,  a*  is 
an  eigenvalue  for  -p  or  equivalently  1-p ,  and  we  have  the  eigenvalue- 
eigenvector  sets 


®  .{♦„(?>} »  p  (2.i) 

®*.  U.*Cy)}’  "p  (2-2) 


Let  us  suppose  further  that  the  undisturbed  flow  has  fore  and  aft  symmetry, 
l.e.  it  is  possible  to  choose  the  origin  so  that 


f (x,y>  -  ¥ (-x,y) 


(2.3) 


Then  Che  perturbation  equation  (1.2)  ia  invariant  under  the 

transformation  t  +  -t ,  x  ♦  -x;  and  we  therefore  have  an  eigenvalue-eigenvector 

aet 


-p  (*•*> 


Taking  the  coaplex  conjugate  of  (2.4)  glvea  the  aet 

^*.{♦*£7)}.  P  (2.5) 

We  conclude  therefore  that  for  undiaturbed  flows  with  fore  and  aft 
symmetry,  and  thla  will  Include  all  the  caaea  considered  in  the  present  paper, 
the  eigenvalues  for  given  p  will  coae  in  pairs  with  a  or  -o*  being  an 
eigenvalue  if  the  other  is.  It  follows  that  the  flow  is  unstable  if  o  is  not 
pure  iaaglnary .  We  also  conclude  that  there  is  a  ayaaetry  about  p  ■  0  or 
p  "  1/2,  since  the  eigenvalue  for  p  la  the  coaplex  conjugate  of  the  one 
for  —p  or  1-p.  Thus  if  o(p)  denotes  the  eigenvalue,  unstable  aodes  satisfy 

Re  a  (p)  ■  Re  o  (-p)  (2.6) 

There  are  also  daaped  nodes  with  growth  rate  -Reo  satisfying  the  same 
equation.  It  must  be  realized,  however,  that  a  nay  be  multivalued  and  the 
equality  in  (2.6)  may  refer  to  different  branches.  It  was  the  failure  to 
realize  this  possibility  that  led  Saffman  and  Schatzman  to  an  incorrect 
assertion  about  the  most  unstable  modes  of  the  Karman  vortex  street . 


When  symmetry  about  the  axis  of  the  motion  exists  it  gives  information  on 
the  parity  of  the  components  of  the  eigenvectors  but  does  not  appear  to  give 
any  further  results  about  the  eigenvalues  themselves.  Reflexional  symmetry 
implies  that  -f (x  +  d,  -y)  describes  the  same  flow  field,  where  d  ■  0  for 
the  single  row  of  symmetrical  double  row,  and  d  ■  1/2  for  the  staggered 
double  row.  The  perturbation  equation  is  invariant  under  the 
transformation  x  x  +  d,  y  ♦  — y ,  -►  -  T ,  and  hence  there  is  an  eigenvalue- 

eigenvector  set 

o,  {e2,rind/%n(-y)},  p  (2.7) 

Por  d  ■  0,  the  eigenvectors  are  therefore  either  even  or  odd  in  y,  and 
for  d  ”  1/2  they  have  a  more  complex  symmetry. 

When  the  vortices  are  small,  the  modes  of  oscillation  will  correspond  to 

those  of  an  array  of  small  circular  vortices,  and  we  expect  therefore  that  the 

disturbances  will  be  of  two  types.  First  will  be  the  cooperative  modes  of 

arrays  of  point  vortices,  for  which  the  characteristic  frequencies  are  of 
2 

order  T It  ,  and  second  will  be  the  oscillations  of  a  single  vortex  for  which 
the  characteristic  frequencies  are  mu/2,  where  u  m  Y  / A  is  the  magnitude  of 
the  vortlclty  inside  the  vortex  of  circulation  T  and  area  A,  and  m  is  the 
angular  wave  number.  These  modes  are  readily  distinguishable,  at  least  for 
small  vortices,  and  provide  a  convenient  way  to  group  the  normal  modes. 

For  the  single  row  (mixing  layer),  the  small  vortex  case  indicates  that 
there  are  two  degrees  of  freedom  for  each  type  of  mode,  and  the  eigenvalues 
for  each  p  will  come  in  pairs,  o^  and  »  where  °2  "  “°i*  the  moc*e  is 
unstable,  and  o^  and  o j  are  pure  imaginary  but  unequal  if  the  mode  is 
stable.  The  symmetry  about  p  *  1/2  also  gives  a  further  two 


eigenvalues  a ,,  *  and  o 2*  for  1-p.  If  unstable,  there  are  then  two  distinct 
branches,  one  with  positive  and  the  other  with  negative  real  part,  which  are 
not  equal  at  p  ■  1/2  and  analytlcity  in  p  and  (2.6)  then  iaplies  that  the 
growth  rate  of  unstable  modes  Is  such  that 

dRe  a 

-  0  for  p  -  1/2  (2.8) 

and  that  |Reo  |  Is  a  maximum  for  the  pairing  mode.  (It  could  In  principle  be 
a  minimum  but  continuous  dependence  upon  the  size  rules  this  out  fo  sell 

vortices,  since  It  Is  a  maximum  for  point  vortices.)  Also  Im(o)  a*  be  zero 

at  p  ■  1/2  for  the  unstable  modes  as  otherwise  the  eigenvalue  wl  r<*t  be 
continuous.  Note  that  p  ■  0  is  a  singular  limit  as  the  disturbance 
wavelength  is  going  to  Infinity,  and  hence  (2.8)  need  not  hold  there. 

In  contrast,  for  the  double  row  (wake),  the  small  vortex  analysis 
Indicates  that  there  are  four  degrees  of  freedom  for  each  mode,  two  for  each 
row,  and  the  eigenvalues  for  each  p  will  come  In  quartets 

<*2»  ®3*  »  wh«r®  if  unstable  o£  ■  and/or  0^  ■  -Oj*. 

There  are  now  four  branches.  The  symmetry  about  p  ■  0  or  p  •  1/2  now 
Implies  that  Oj*(p)  •Oj(-p),  etc.  (o^*  (p)  -02  ("P)  violates  the 
continuity  of  Im(o)  at  p  ■  1/2  unless  the  Imaginary  part  vanishes,  and  this 
Is  contradicted  by  the  small  vortex  results.)  It  does  not  follow  therefore 
that  Re  o^  Is  symmetrical  about  p  ■  0  or  p  ■  1/2  when  there  Is 
Instability.  For  p  ■  1/2,  the  quartet  of  eigenvalues  will  be  such 
that  0,  0*,  -0 ,  and  -0*  are  all  eigenvalues,  but  the  real  and  imaginary 
parts  that  coincide  come  from  different  branches,  and  (2.8)  does  not 
necessarily  hold.  Contrary  to  the  statement  by  Saffman  and  Schatzman  and  in 
agreement  with  Rida's  results,  the  p  ■  V2  disturbance  will  not  be  the  most 


unstable  disturbance. 

In  figure  1,  we  show  a  sketch  of  the  dependence  of  the  four  complex 

growth  rates  on  p  for  the  cooperative  modes  for  zero  area  and  general  ic  .  In 

this  special  case,  the  system  Is  degenerate  and  there  Is  symmetry  about 

both  a  ■  0  and  p  ■  1/2.  In  addition,  for  the  exceptional 

value  h/i  ■  k  -  k  ,  the  circle  of  Instability  reduces  to  a  point  and  the 
c 

imaginary  values  of  a  meet  in  a  cross.  Figures  2,  3  and  4  show  how  these 

curves  may  change  when  the  area  of  the  vortices  is  finite.  Figure  2  is  an 

example  where  the  street  remains  unstable  for  p  -  1/2,  but  the  pairing 

Instability  Is  not  the  most  unstable  disturbance.  Figure  3  Is  a  sketch  of  a 

possible  behavior  in  which  the  street  is  unstable,  but  the  pairing  instability 

with  p  -  1/2  is  stabilized  by  finite  area.  Figure  4  is  the  case  where 

finite  area  has  stabilized  the  street.  It  can  be  seen  that  the  symmetry  of 

reflection  about  a  ■  0  and  p  ■  1/2  is  satisfied,  but  the  symmetry  about  the 

axes,  which  is  true  for  zero  area,  breaks  down  when  the  area  is  finite. 

For  k  away  from  tcc>  it  is  clear  that  figure  1  will  first  change  into 

figure  2  when  the  area  becomes  finite  and  may  later  change  into  figures  3  and 

4  when  the  area  becomes  larger.  Note  that  these  results  depend  upon  fore  and 

aft  symmetry,  but  the  sizes  of  the  vortices  in  one  row  need  not  be  the  same  as 

those  in  the  other  row.  The  task  is  to  determine  whether  the  special  case  of 

figure  1  for  ic  ■  changes  into  figure  2,  3  or  4  when  a  *  A/A  is  finite  for 

values  of  <  very  close  to  k  .  There  is  also  a  fourth  possibility  that  for 

c 

finite  area  there  is  an  exceptional  value  of  k,  which  should  depend  upon  area, 
such  that  figure  3  is  degenerate,  i.e.  the  ovals  shrink  into  points  and  the 
imaginary  values  meet  in  crosses. 

The  high  frequency  shape  modes,  which  do  not  exist  for  point  vortices, 
are  found  to  be  stable,  provided  the  vortices  are  not  too  large,  and  their 
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Figure  1. 

Sketch  of  dependence  of  real  and  Imaginary  parts  of  the  growth  rate  on  wavenumber 

p  for  zero  area  and  a  typical  value  of  k  •  For  ic  ■  k  ,  the  oval  shrinks  to  a  point 

c 

and  the  imaginary  parts  form  two  crosses. 


Figure  2 . 

Sketch  of  dependence  on  p  of  an  unstable  quartet  of  eigenvalues  for  fixed  area  and 
aspect  ratio.  The  pairing  mode  is  unstable,  but  is  not  the  most  unstable 
disturbance.  Note  the  form  of  the  symmetry. 


Figure  3. 

Sketch  of  dependence  on  p  of  an  unstable  quartet  of  eigenvalues  when  the  pairing 
mode  Is  stable*  A  degenerate  case  may  exist  in  which  the  ovals  shrink  to  two  point 
and  laaglnary  parts  fora  two  crosses. 


Figure  4. 
Sketch  of  de 


cr  on 


Wh,n  *  «.n.  occur. 
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oscillation  frsquaoclas,  which  will  bs  modified  by  tha  Interactions,  are  as  In 
figure  4. 

We  now  concentrate  on  the  staggered  double  row  and  discuss  the  detailed 
problea  of  what  should  be  calculated  and  how  we  shall  proceed. 

3.  The  Karaan  street  of  aaall  vortices 

The  coaplex  growth  rates  a  will  ba  the  eigenvalues  of  an  Infinite  systea 

M*  -  04  (3.1) 

where  M  will  have  an  expansion  of  the  fora 

^4M-1  +M0+j4M1  +frM2  +,,,1  (3-2) 

Invariance  with  respect  to  the  sign  of  1  requires  that  the  expansion  of  N 
goes  In  powers  of  a  and  not  In  powers  of  /a  .although  the  detailed  solution  Is 
an  expansion  in  1 H  .  The  aatrlces  depend  upon  the  aspect  ratio  k  and  the 

subharaonlc  wavenuaber  p.  The  aatrlx  M_j  which  produces  shape  deforaatlon 
nodes  Is  given  by  the  slaple  calculation  for  the  oscillations  of  an  isolated 
circular  vortex,  while  Mq  which  produces  the  cooperative  oscillations  of  the 
centroids  Is  found  froa  the  Karaan  calculation  of  point  vortices.  M^,  etc., 
have  to  be  found  by  laborious  algebra.  The  eigenvectors  of  M_}  and  Mg  are 
Independent.  In  the  absence  of  degeneracy,  which  occurs  when  two  or  more  of 
the  eigenvalues  are  equal,  the  eigenvalues  will  be  analytic  functions  of  the 
paraaeters. 

Proa  the  general  aymaetry  arguments  of  section  2,  we  know  that  the 
eigenvalues  will  appear  in  quartets,  which  will  satisfy  quadratic  equations 


(3.3) 


r  2  rz 

[fl  -i-rSt(p,K,oi)]  - -T  H.  (p  ,ie  ,a  ) 
i  t 

2 

[o  -  i  -y  S3(p,ie  ,o)]2  -  H3(p,tc  ,o)  (3.4) 


where  S^,  83,  Hj  end  H3  are  real  and  dimensionless.  The  symmetry  about 
p  ■  0  or  p  ■  1/2  Implies 

S3(p,*)  -  -S1(-p,*),  H3(p>*)  -  31(-p,*)  (3.5) 


Non-trlvial  degeneracy  occurs  If  and  when  and  H3  vanish,  and 

changes  of  stability,  where  o  changes  from  complex  to  pure  Imaginary  can  only 
occur  at  such  points.  With  regard,  therefore,  to  the  main  problem  of  the 
stabilization  of  the  Karman  vortex  street  by  finite  area  vortices,  the 
Interest  Is  In  the  functions  and  H3  and  their  zeros. 

These  functions  are  known  In  the  limit  a  *  0,  which  can  be  thought  of  as 
either  finite  size  vortices  Infinitely  far  apart,  as  Is  appropriate  for  the 
shape  modes,  or  point  vortices  separated  by  a  finite  distance,  as  Is 
appropriate  for  the  cooperative  modes.  It  appears  from  our  present  numerical 
calculations  and  those  of  Saffman  and  Schatzman  that  and  H3  are 

negative  for  the  shape  modes  when  a  >  0  (but  not  too  large)  which  therefore 
remain  stable,  and  we  leave  these  modes  for  the  present  and  concentrate  on  the 
cooperative  modes.  In  this  case, 

2  2  2 

^(p.x.O)  -  H3(p,K,0)  ■  -j—  [  2p(l-p)  -  sech  (me)] 

2  2  4 

-  t-[ cosh(2mep)  -  2p  cosh(«K)  cosh (me (l-2p) )  ]  /cosh  (me) 

(3.6) 


The  equality  of  and  H3  when  a  •  0  means  that  for  point  vortices  the 


eigenvalues  occur  in  pairs  such  that  -o  is  an  eigenvalue  if  a  is.  This  is 

also  a  consequence  of  the  fact  that  point  vortices  fora  a  Hamiltonian  system 

with  x  and  y  coordinates  as  conjugate  variables.  It  is  Important  to 

realise  that  this  additional  symmetry  is  a  consequence  of  a  zero  area 

degeneracy  and  is  broken  at  finite  area. 

The  pic  plane  for  a  ■  0  is  divided  up  into  regions  of  stability  and 

instability  and  the  boundaries  look  like  a  diagonal  cross  with  center  at 

p  ■  1/2  and  k  ■  k  ,  where  cosh(«K  )  •  /T  (see  Saffman  and  Schatzman  (1982) 
c  c 

figure  2).  For  k  not  close  to  kc,  the  effect  of  finite  area  is  to  displace 
the  boundaries  of  the  stability  regions ,  and  the  results  are  not  of  pressing 
Importance.  For  sswll  area,  the  main  interest  lies  in  the  behavior  of  p 
close  to  1/2  and  <  nearly  equal  to  k£.  Let 

p'  -  p  -  1/2,  k*  -  x  -  *c  (3.7) 

then  we  need  the  expansion  of  and  H3  for  small  p',  k'  and  a. 

The  numerical  results  of  Saffman  and  Schatzman  for  p'  ■  0,  our  present 
arguments  and  results  (see  last  paragraph  of  §4)  and  the  analysis  of  Klda 
indicate  that  there  are  no  terms  in  the  expansion  linear  in  a,  and  that  the 
dominant  terms  of  the  expansion  are 

Hl/3(p,K,0)  -  c1ic,2+  c2p,2+  CjO2*  ’+/-  c2u2p'+  c^a4+  •••  (3.8) 

The  first  coefficient,  c^,  is  positive  and  the  second  one,  C2,  is 
negative.  Both  are  given  analytically  by  the  point  vortex  theory.  The  third 
and  fourth  are  given  analytically  by  the  Kida  expansion.  The  last 
coefficient,  C5,  can  be  deduced  from  the  Saffman  and  Schatzman  calculation, 


which  demonstrated  that  it  waa  negative  (and  aleo  givea  C3) •  Theae  numerical 

3  3  3 

reaults  exclude  ana  tern.  Ana  p'  or  a  k'  tern  could  be  present,  but  la 

unimportant  conpared  with  the  a^  tern  In  ita  effect  on  the  atabllity 

4 

boundary.  Notice  that  Rida's  expanaion  doea  not  Include  all  a  terne 
consistently,  alnce  they  are  not  Included  In  the  dynamics,  and  hence  hla 
expreaalon  for  the  fifth  tern  la  Incomplete.  Hla  quantitative  conclusions 
about  the  stabilisation  near  the  critical  point  are  therefore  Invalid. 

The  numerical  values  of  theae  coefficients  determined  by  the  earlier  work 
are  as  follows.  From  the  Rarman  expressions, 

cx  -  12.1761,  c2  -  -0.700541  (3.9) 

From  the  results  for  the  stability  boundary  for  p  -  1/2  given  in  equation 
(3.19)  of  the  Saffman  and  Schatxman  paper,  which  relate  the  values 
of  a  and  ic  for  which  Hj  or  H3  vanish  when  p'  •  0,  we  have  the 

estimates 

I 

I 

I 

cS3  -  -12.88,  cS5  -  -11.60  (3.10) 

i 

I 

|  From  the  analytical  expression  given  by  Rida's  equation  (3.29),  and  the 

formulae  of  his  (3 .10)-(3 .11) ,  we  find  that 

cK3  -  -13.52,  c\  -  6.49,  cK5  -  -12.33  (3.11) 

It  can  be  seen  that  the  C3  and  C5  are  in  reasonable  agreement.  This  is 
consistent  with  Rida's  statement  that  hla  reaults  agree  with  those  of  Saffman 
and  Schatzoan  for  p  »  1/2,  as  will  be  discussed  in  more  detail  at  the  end  of 
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th«  section,  see  equations  (3.20)-(3.22) . 

The  stabilisation  of  the  street  Is  determined  for  small  area  by  (3.8). 
We  complete  the  squares  and  writs  the  right-hand-side  of  (3.8)  as 

CiCk'+c^2)2  +  c2(p'+/-  c4«2)2 
"**1  ""**2 

(3.12) 

+  o4(c5-c23/4Cl-c24/4c2) 

The  stabilisation  or  destabilisation  of  the  street  by  finite  area  therefore 
depends  upon  the  sign  of 


D  -  c5  -  c23/4Cl  -  c24/4c2  (3.13) 

If  D  Is  positive  the  street  Is  made  unstable  In  the  vicinity  of  the  critical 

point  by  finite  area.  If  D  Is  negative  the  street  Is  stabilised  for  a  range 
2 

of  k  of  order  a  by  finite  area. 

If  D  ■  0,  then  the  critical  point  remains  an  isolated  value  of 

stability  but  the  aspect  ratio  of  the  stable  street  Is  eltered  by  a  value  of 
2 

order  a  •  In  this  case.  It  would  be  necessary  to  go  to  yet  higher  order 
In  a  to  determine  the  stability  question. 

For  a  fixed  a,  the  stability  boundaries  in  the  ic'-p'  plane  are  given  by 
the  vanishing  of  (3.12)  which  gives  two  hyperbolae.  The  interiors  of  the 
hyperbolae,  marked  by  cross-hatching,  contain  the  unstable  modes.  There  are 
two  main  possibilities  depending  upon  the  values  of  the  c  coefficients. 
Figure  5  shows  the  stability  boundaries  for  p  close  to  1/2 
and  k  near  «c  when  D  Is  negative.  There  is  a  region  of  stability  bounded 
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Figure  5. 

Stability  boundary  in  the  k'  -  p'  plane  when  finite  area  stabilizes.  Hyperbolae 
would  cross  at  centers  if  D  ■  0.  Cross  hatched  areas  denote  Instability  regions 
for  modes  1  and  3.  Numbers  in  parentheses  show  corresponding  figure  for  growth  rate 
as  function  of  p. 


2  2 

by  the  two  hyperbola*  with  centers  at  k*  ■  Xa  and  p*  *  *  ua  ,  where 


X  »  -c« 


W  -  -c. 


(3.14) 


The  distance  between  the  vertices  giving  the  range  of  k  for  which  there 


Is  stability  la 


fiic  -  (-D/c1)1/2a2 


(3.15) 


In  this  case,  finite  area  stabilises  the  street  for  a  rang*  of  aspect  ratio. 


On  the  other  hand,  If  0  Is  positive,  then  the  stability  boundarl« 


as  shown  In  figure  6.  Now,  there  Is  always  a  value  of  p  which  give*  an 


unstable  aode  and  finite  area  destabilizes  the  street  by  removing  the  neutral 


case .  If 


M  >  (-D/c2)J 


(3.16) 


there  Is  a  lens  of  stability  along  the  p'  ■  0  axis  and  the  boundaries  are  as 


shown  In  figure  6.  The  pairing  aode  is  then  stabilized,  but  there  Is  a  range 


of  unstable  values  of  p  of  width 


.  ,  .1/2  2 
fip  ■  (-0/c2)  a 


(3.17) 


If  (3.16)  Is  not  satisfied,  the  hyperbolae  overlap  and  there  is  no 


stabilization  of  the  pairing  mode. 


If  D  ■  0,  the  hyperbolae  degenerate  into  crosses  meeting  at  the  points 


given  by  (3.14)  and  the  street  Is  stable  for  one  special  value  of  k  equal 


V.-.  .•>. 


--.‘A  •  1 
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Figure  6. 

Stability  boundary  In  the  k'  -  p'  plane  when  finite  area  destabilizes.  In  case 
shown,  the  pairing  instability  p '  •  0  is  stabilized. 


2 

to  <  +  Xo  . 

c 

Using  Rida’s  values  given  by  (3.11)  we  find  that  D  ■  **1.04.  Using 

Saffman  and  Schatzman's  estimate  of  03  and  03,  and  Rida's  value  of  04, 

we  have  D  ■  0.037.  There  Is  no  reason  why  these  estimates  should  not  be  of 

opposite  sign,  because  as  explained  above  there  Is  an  error  in  Rida's  cj*  It 

appears,  therefore  that  0  is  small  and  the  question  of  stabilization  by 

finite  area  will  be  delicate.  Either  a  consistent  perturbation  expansion  of 
4 

order  a  should  be  carried  out,  or  a  further  numerical  Investigation  be 
attempted.  The  algebra  entailed  in  the  former  is  horrific  (Saffman  and 
Schatzman  did  consider  carrying  out  such  an  expansion  for  the  simpler 
p  ■  1/2  case,  but  decided  In  favor  of  a  numerical  calculation)',  and  it  was 
again  decided  to  go  to  a  numerical  approach. 

We  discuss  in  the  next  two  sections  our  method  for  calculating  the  steady 
shapes  and  stability  boundaries,  which  we  use  to  investigate  the  stability  of 
the  cooperative  modes.  It  also  gives  the  high  frequency  shape  modes,  and  has 
the  advantage  that  it  Is  easily  applied  to  a  wider  class  of  problems.  With 
respect  to  the  street  stability,  the  numerical  results  Indicate  that  D  Is 
actually  zero  to  within  the  limit  of  numerical  accuracy.  Rida's  results  that 
there  Is  a  range  of  k  for  which  finite  area  stabilizes  disturbances  of 
arbitrary  p  appears  therefore  to  be  Incorrect. 

On  the  other  hand,  Rida's  expansion  should  give  the  coefficient 

2 

of  a  exactly  in  closed  form  for  all  p  and  k.  Away  from  the  critical  point, 

there  will  be  values  p  ,  k  at  which  the  stability  of  point  vortices 

8  8 

changes.  To  determine  the  effects  of  finite  area  for  these  values,  we  require 
an  expansion 

Uj^P.K.a)  •  c1(k-ks)  +  c2(p-pg)  +  c3<x2  (3.18) 
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and  a  similar  expansion  for  H3.  Here,  c ,  and  C2  ere  known  functions  of 

Pa  which  vanish  for  pg  ■  1/2,  and  the  same  is  true  for  C3  which  is  given 

by  Kida's  results*  For  ps  *  1/2,  (3.18)  is  sufficient  for  the  determination 

of  the  stability  properties  for  small  a>  However,  the  coefficients  in  (3.18) 

vanish  when  pfl  *  1/2,  and  there  the  expansion  (3.8)  is  required. 

The  discussion  of  this  section  shows  that  from  Rida’s  expressions  for 

2 

H}  and  H3  correct  to  0(a  )  combined  with  the  numerical  results  which  give 

the  coefficient  C5,  we  can  construct  an  expression  for  the  effect  of  finite 

2 

size  on  the  stability  boundaries  which  is  uniformly  valid  to  0(a  )  for  all 
p  and  <  .  We  have  for  the  cooperative  modes 

Hi(p,<,a)  ■  "linear  theory" 

+  "Rida"  a2  (3.19) 

+  "numerical" 

with  H3  given  by  (3.5).  Rida's  analysis  also  gives  Sj  and  S3  correct 
to  0(a2) . 

For  the  special  case  of  disturbances  with  p  *  1/2,  p'  -  0,  Saffman  and 
Schatzman  found  that  for  small  area  the  street  was  stabilized  for  aspect 
ratios  in  the  range 


-0.583  a2  <  <’  <  1.644  a2  (3.20) 

According  to  (3.8),  the  range  of  values  is 

-b.a2  <  k  '  <  b0a 


(3.21) 


C3  *  ci(bi  -  b2>,  C5  ■  -cjb^ 


(3.22) 


Klda's  values  of  03  and  05  give  b^  ■  0.60,  b2  "  1.70;  which  agree 
reasonably  with  the  values  found  by  Saffman  and  Schatzman. 

4.  Calculation  of  steady  shapes 

We  describe  the  method  here  In  the  context  of  the  staggered  double  row. 
Modifications  for  other  flow  configurations  are  simple  and  easy  to 
Implement.  Saffman  and  Schatzman  employed  a  boundary  Integral  method.  This 
has  the  advantage  that  It  can  be  employed  for  finite  area,  but  It  suffers  from 
the  disadvantage  of  being  rather  hard  to  Implement  for  disturbances  of 
arbitrary  wavenumber  (although  this  possibility  is  under  investigation)  and 
also  being  hard  to  check  analytically.  The  method  we  employ  here  follows  one 
suggested  by  Dr.  Javier  Jimenez  and  Is  based  on  an  expansion  in  inverse  powers 
of  1  .  This  can  be  thought  of  as  an  expansion  In  core  radius  or  /a ,  but  we 
choose  to  keep  the  area  of  the  vortices  finite  and  suppose  the  separation  t 
Is  large,  with  h/1  kept  fixed.  In  principle,  the  algebra  can  be  done  and  a 
consistent  perturbation  expansion  developed,  but  this  is  beyond  our  present 
resources  (even  with  the  aid  of  symbolic  manipulation  programs)  and  we  do  the 
algebra  arithmetically.  Unfortunately,  the  detailed  approach  does  not  allow 
the  retention  of  high  order  terms  In  a  completely  uniform  manner,  so  checking 
with  Kids' s  perturbation  results  suffers  from  uncertainty. 

The  shape  of  a  vortex  in  the  first  row  is  assumed  to  be  described  by  the 
exterior  conformal  map  of  the  unit  circle  |c|  ■  1; 

z  -  RjC{l  +  a^/z  +  a2/c2  +  •••}  (4.1) 


■j'j 
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Vortices  in  the  second  row  have  the  expansion 


z  -  R2c{1  +  b1/c  +  b2/?2+  •••}  (4.2) 

These  series  may  fail  to  converge  before  the  solutiona  actually  break  down  by 
vortices  coming  into  contact,  but  we  assume  that  there  is  convergence  for  a 
finite  range  of  area. 

Our  first  task  is  to  calculate  the  velocity  field  Induced  by  a  single 
vortex  containing  vorticity  of  shape  (4.1).  To  do  this,  we  employ 
(following  Jimenez)  the  so-called  Schwarz  functions,  which  are  analytic 
functions  of  z  equal  to  z*  on  a  contour.  To  find  the  Schwarz  function  for 
the  contour  (4.1),  we  expand  the  expression  conjugate  to  (4.1),  remembering 
that  5*  ■  1/5  on  the  contour, 


z*  -  R*1/5{1  +  a*tc  +  a*2c2+  •••}  (4.3) 

as  a  Laurent  series  in  z,  giving  an  expansion 

**  ■  8q/z  +  8i/z2  +  82/z^+,*,+fo  +  *1*  +  f2z^+***  (4.4) 


The  complex  velocity  u  -  iv  induced  by  the  vortex  is  analytic  outside  the 
vortex,  and  such  that  u  -  iv  +  lu^z*/2  is  analytic  inside  the  vortex.  It 
follows  by  inspection  that  the  exterior  velocity  field  is 


u  -  iv 


-iw 


!  80  gi  S2 

-  ll  +  ~T  +  T  + 


(4.5) 


The  interior  velocity  field  which  can  be  expressed  in  terms  of  the  f_, 


etc . , 


Is  not  required.  It  remains  to  describe  how  the  gn  are  found  In  terms  of 
the  an. 

To  do  this,  ve  multiply  (4.4)  by  sn  and  Integrate  around  a  contour  In 
the  c-plane,  giving 

8n "  sbr*  **  *“  **  "  SFT*  **  *°  If* 

■  coefft  of  In  product  z *  *n  ^  (4.6) 

when  (4.1)  Is  substituted  for  z  and  (4.3)  for  s*. 

Similarly. 

g*n  ■  -coefft  of  In  product  z  z*a  — V  (4.7) 

the  minus  sign  arising  from  reversal  of  the  sense  of  the  contour. 

There  are  alternative  expreslons  for  the  gn  In  terms  of  moments  of  the 
shape  of  the  vortex.  Since 

U  ”  lv  "  F“l  //  ~^T$r  (4.8) 

It  follows  that 

v  gn-JJ  *'°dx’dy’  (4.9) 

where  the  integral  is  over  the  vortex.  In  particular,  gg  is 

the  area /» ,  and  g^  is  proportional  to  the  displacement  of  the  centroid 

relative  to  the  z  origin. 

Similarly,  we  obtain  an  expression  for  velocity  Induced  by  the  vortex 


vv 


u  -  iv 


_ *  /  «  4  f. 


(4.10) 


where  the  h.  ere  related  to  the  b_  by  the  fornula  enelogoue  to  (4.6) 


We  alao  obtain  the  conjugate  velocity 


+  Iv 


*“l  /8*°  x  8*1  .  8*2 

~r 

**  ** 


(4.11) 


and  similarly  for  the  second  vortex,  where  the  g*a  are  related  to  the  an 
and  a*n  by  (4.7). 

An  expression  can  now  be  written  down  for  the  complex  velocity  produced 
by  the  two  staggered  rows.  The  origins  in  the  first  row  will  be  supposed  to 
be  at  the  points  ml ,  <  m  <  • ,  and  In  the  second  row  at  the 

points  ml  -  A(d  +  lx).  The  separation  1  need  not  be  real;  Its  phase 
determines  the  angle  between  the  x-axis  and  the  direction  of  the  rows. 

Indeed,  It  proved  convenient  in  the  actual  calculations  to  take  t  to  be  pure 
imaginary.  The  relative  stagger  d  and  the  aspect  ratio  x  are  real.  For  the 
Karman  vortex  street,  d  ■  1/2.  For  the  symmetrical  double  row,  d  -  0. 
Arbitrary  values  of  d  may  be  considered,  but  we  shall  not  do  so.  Then  at  a 
point  with  coordinate  z  relative  to  a  vortex  in  the  first  row 


-ia).  -  g„ 

•-"•-l 1  I  l  -f->i 

n-0  in--®  (z-o£  ) 


-i/Do  •  « 

z  r  r 


l  i  — 2 - 

n-0  m— •  (z-mi  +!(d  +  lx))1 


+  U  -  IV, 


(4.12) 


and  at  a  point  with  coordinate  z  relative  to  a  vortex  In  the  second  row 


U  -  iv  -  -5—  l  £ 


n«0  n— •  (*-*t ) 


• 

r  I  l 

1  n-0  ■- 


—  (x-m4-i(d+iic)) 


+  0  -  iV.  (4.13) 


U  and  V  are  components  of  "ho  velocity  at  infinity  which  balances  tha  aalf 
induced  notion  of  tha  atraat  and  brings  the  vortices  to  rest. 

Expressions  for  the  conjugate  velocities  u  +  iv  are  obtained  by 
replacing  i  with  -i ,  gQ  with  g*n,  hn  with  h*n,  and  i  with  1*. 

For  steady  notion,  the  shapes  of  the  vortices  are  given  by  the  condition 


that  their  boundarl* 


itreaallnea,  l.e.  u  -  iv  is  parallel  to  the 


tangent  ds/ds.  This  is  expressed  by  the  equation; 


«  ft)  *  (u-iv)c  45.+  (u+iv)c  -il*-  0  (4.14) 


when  |c|  ■  1 ,  to  be  satisfied  on  the  first  and  second  vortices,  l.e.  with 
u-  iv  given  by  (4.12)  and  s  given  by  (4.1),  and  u  -  iv  given  by  (4.13) 
and  z  given  by  (4.2),  and  the  appropriate  conjugate  expressions, 
respectively.  This  provides  equations  for  the  unknown  coefficients  Rj_,  an, 
r2»  bn*  translation  velocity  Q  ■  D  -  IV  and  Q*  ■  0  +  IV. 

It  is  expected  that  locally  unique  (Isolated)  solutions  will  exist  if  the 
areas  of  the  vortices,  the  positions  of  the  centroids,  and  the  strengths  of 
the  vortices  are  given.  The  question  to  be  discussed  now  is  the  actual 
procedure  to  be  employed. 

The  straightforward  way,  and  the  one  that  is  implicit  in  the  present 


formulation,  is  to  develop  a  series  expansion  in  a,  the  solution 


for  a  •  0  being  known.  It  Is  convenient  to  suppose  that  the  vortices  have 
finite  areas  and  strengths,  and  that  the  Halt  a  ■  0  is  obtained  by 
letting  i  ♦  •  .  Then  we  can  develop  the  velocity  by  expanding  the  teros  in 
(4.12)  and  (4.13)  and  their  conjugates  as  power  series  in  1/1  and  calculate 


the  unknowns 


rles  in  the  same  variable.  This  is  feasible  in  principle. 


but  iapos8ible  in  practice  at  present,  especially  when  it  is  remembered  that 

->9 

the  stabilization  problem  requires  retention  of  terms  of  order  1  in  order  to 

A 

retain  consistently  effects  of  order  a  .  An  alternative  approach,  and  the  one 
that  we  employed,  is  to  use  the  computer  to  do  the  algebra  arithmetically  as 
follows . 

We  pick  two  Integers  N  and  L.  The  series  (4.1),  (4.2)  and  their 
conjugates  are  truncated  to  keep  the  first  N  +  1  terms,  l.e.  we  Include  aN, 
bjj,  a*N  and  b*N.  Newton  Iteration  will  be  employed  to  find  the  values  of 
these  (in  general)  complex  numbers.  For  this  purpose  we  then  calculate 
numerically,  starting  with  a  first  guess  for  the  2N  +  2  unknowns  Rlt 
ai,..., aN,  R*i,  a*i,...,a*u,  from  (4.6)  and  (4.7),  the  values  of  gQ, 
8l».*.,gfl,  g*o«  g*i » •  •  •  »g*jj •  In  fact,  go  and  g*o  are  put  equal  to  the 
given  area  of  the  vortex  divided  by  it ,  and  the  difference  between  the 

values  calculated  from  (4.6)  and  (4.7)  are  used  as  a  check  on  the  accuracy  of 
the  calculation.  Similarly,  we  calculate  the  hQ  and  h*„.  For  reasons  that 
are  important  for  the  stability  calculation  and  will  be  made  clear  later,  we 
treat  quantities  and  their  conjugates  as  independent  variables,  e.g.  Rj 
and  R*^  are  regarded  as  Independent  complex  numbers.  Of  course,  the 
solution  will  not  be  physically  sensible  unless  the  final  answers  are  complex 
conjugates  of  each  other.  So  far,  we  then  have  4N  +  4  complex  unknowns. 

We  now  expand  (4.12)  as  a  power  series  in  1/1  retaining  terms  of 
order  1  **.  This  gives  a  triple  series 


»v-v. 

a.  -  a’  .  •«  .  *  -  -  . 

_•  '  t 


i 
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V  * 


where  B(p,n)  la  tha  eoafflelaot  of  zp  In  tha  binoalal  expansion  of 
(1+*)”°.  Tha  vortlcitles  ara  ta.  -  T/A^  and  0*2  “  -T/Aj  •  Tha  auaa  with 
raapaet  to  a  can  ba  dona  in  cloaad  fora,  aea  appendix  B.  Ha  now  aubetitute 
eha  truncated  axpanalon  (4.1)  into  (4.15)  and  obtain  u  -  lv  evaluated  on  tha 
boundary  of  tba  vortex  aa  a  aarlaa  of  poaltlva  and  negative  powers  of  c  • 
Differentiation  of  (4.1),  ate.,  glvea  a  aarlaa  in  C  for  dz/d?  ,  etc.,  and 
aubatltution  into  (4.14)  glvea  a  aarlaa  of  poaltlva  and  negative  powers 
of  C  with  coaplax  coefficients  which  are  analytical  functions  of  d 
and  k,  contain  inverse  powers  of  t,  and  are  nuaerlcal  functions  of  the 
4H  +  4  unknowns,  Rj,  an,  etc.  They  also  contain  the  two  unknown  coaplex 
velocities  Q  ■  U  -  IV  and  Q*  •  0  +  IV,  which  will  also  be  regarded  aa 
independent  coaplex  unknowns. 

To  obtain  equations  for  these  4N  +  6  unknowns,  we  calculate  the 
coefficients  of  ?n  In  the  expression  (4.14)  for  values  of  n  in  the 
range  -N  <  n  <  N,  for  vortices  in  the  first  and  second  row,  and  denote  these 
4N  +  2  coaplex  quantities  (2N  +  1  for  each  row)  by  E}(n)  and  E2(n).  He 
obtain  4n  +  2  equations  for  the  unknowns  by  requiring 

E1(n)  -  0,  E2(n)  -  0,  -N  <  n  <  N  (4.16a,b) 

He  have  not  yet  said  anything  about  the  position  of  the  origins  inside 


the  vortices,  from  which  the  seperetlon  k  Is  measured.  It  Is  convenient  to 
require  that  the  centroids  coincide  with  the  origins.  These  constraints  upon 
the  variables  are  equivalent  to  the  four  equations 

gl  “  g*i  ■  hi  •  h*i  ®  0  (4.17) 

Two  further  equations  are  found  from  the  condition  that  the  vortices  have 
given  area.  Defining  the  area  by 

A1  "  *  RlR*Jl  “  *2**2  “  2*3**3  “  •••“(N-l)aMa*K}  C**1®) 

and  a  similar  expression  for  A2»  which  are  correct  when  starred  quantities 
are  complex  conjugates,  and  which  are  Identical  to  gQ  -  g*g  “  fa  in  this 
case,  we  have  two  further  equations. 

Thus  altogether,  we  have  4N  +  8  equations  for  4N  +  6  unknowns.  The 
equations  are,  however,  not  Independent.  Firstly,  E^(0)  and  E2(0)  are 
Identically  zero,  and  secondly  there  are  two  independent  constraints  on  the 
E^Cn)  and  E2(n).  The  reasons  are  as  follows. 

On  the  surface  of  the  vortex 

(u  -  iv)dz  -  (ug  +  iun)ds  (4.19) 

where  u8  and  un  are  the  tangential  and  normal  components  and  ds  Is  the 

element  of  length.  Also,  ?  •  exp(18).  It  follows  that  2u  ds  ■  R  (c )  d8  . 

n 

Now  the  velocity  field  that  we  construct  above  is  an  exact  solution  of  the 
Euler  equations  and  hence  must  automatically  satisfy  the  equation  of 
continuity.  Thus,  automatically 


£uQds  -  0  •  Jq  R  (Ode  -  2vE(0) 


(4.20) 


and  hence  the  vanishing  of  Ej(0)  and  E2<0)  are  Juat  consistency  checks  and 
not  Independent  equations. 

Neat,  consider  the  consistency  condition  (A6)  derived  in  Appendix  A. 
Since, 


/  /  (u  A  iv)dxdy  ■  P(x  ±  iy)  uQds 
we  can  rewrite  this  equation  for  constant  vorticlty  as 

w  *  z*)  R  (?  )dd  +u>2f2(x  ±  **)  R  (C)d0  -  0 

(4.21) 

Thus  we  have  two  constraints  relating  the  components  of  R  (?)  on  the  two 
vortices  which  are  imposed  by  klneaatica.  Since  z  is  given  by  (4.1),  etc., 
the  relations  are  between  all  the  R|(n)  and  E2(n),  but  since  the  dominant 
terns  in  the  expansions  are  the  leading  ones,  it  is  expected  that  it  suffices 
to  throw  away  the  equations  for  E2(l)  and  E2(-l).  The  existence  of  (4.21) 
will  ensure  that  the  converged  solution  will  make  these  quantities  zero 
automatically ,  and  in  fact  the  accuracy  of  the  solution  can  be  monitored  by 
their  values. 

To  summarise,  we  now  have  4N  +  4  independent  complex  equations  from 
(4.17)  (4  equations),  (4.18)  (2  equations),  (4.16a)  with  n  ■  0  excluded  (2N 
equations),  and  (4.16b)  with  n  ■  0,  n  ■  1  and  n  *  -1  excluded  (2N  -  2 
equations).  However,  we  still  have  4N  +  6  complex  unknowns  R},  anj,  R2» 


a  2>  their  conjugates,  and  the  two  complex  velocities.  Two  further  equations 


are  therefore  required*  These  ere  obtained  by  noting  that  the  phase  of  c  is 
arbitrary,  and  that  this  can  be  fixed  by  specifying  the  phases  of  and 

Rj*  He  do  this  by  requiring  that  in  the  final  solution  R}  and  R2  should 
be  real,  which  is  imposed  by  the  two  complex  equations 

Rx  -  R*l,  R2  -  R*2  (4.22) 

The  number  of  equations  and  unknowns  is  now  equal.  The  test  that  the 
counting  has  been  done  consistently  is  the  non-vanishing  of  the  Jacobian  of 
the  system.  (The  constraints  embodied  in  (4.21)  were  in  fact  discovered 
numerically  by  the  vanishing  of  the  Jacobian  when  the  n  -  +/-1  equations  of 
(4.16b)  were  not  excluded.) 

For  given  values  of  the  parameters,  a  first  guess  was  taken,  usually  the 

circular  vortex  approximation,  and  solutions  were  calculated  by  solving  the 

system  with  Newton  iteration.  It  is  necessary  to  decide  on  suitable  values 

of  N  and  L.  Taking  L  ■  1  gives  circular  vortices,  since  the  induced 

velocity  in  the  neighborhood  of  a  vortex  is  then  constant.  To  calculate 

deformation  consistently  to  order  a,  it  is  necessary  to  take  L  -  2  and  it  is 

then  sufficient  to  take  N  ■  2  which  gives  vortices  of  elliptical  shape. 

This  corresponds  exactly  to  what  can  be  called  the  elliptical  vortex 

approximation  (Saffman  and  Szeto  1981).  The  exterior  velocity  Induced  by  a 

2 

uniform  vortex  of  elliptical  shape  differs  by  terms  of  relative  order  a  from 

that  of  a  circular  vortex  of  the  same  circulation,  so  at  this  order  the 

cooperative  behavior,  i.e.  translation  velocity  and  stability  properties,  will 

be  the  same  as  for  point  vortices.  This  explains  why  there  are  no  terms 

linear  in  a  In  (3.8).  The  first  significant  effect  of  area  comes  with  L  - 

2 

3,  which  gives  a  consistent  calculation  of  0(a  )  effects.  The  value  of  N 


should  b«  at  laaat  L,  and  there  la  little  point  in  taking  larger  values . 

The  accuracy  of  the  solution  can  be  tested  by  comparing  vlth  the  results  of 
Saffman  and  Schatzman  (1982)  and  the  steady  solutions  correct  to 
order  a  'given  by  Kids.  tJnfortunately(  the  numerical  method  does  not  Ignore 
smaller  terms,  and  results  obtained  with  L  -  5  cannot  agree  exactly  frith 

4 

Kids' s,  but  tha  difference  should  be  of  order  a  •  For  the  investigation  of 
stability,  at  least  near  the  critical  value  of  ic  »  iec,  it  is  necessary  to 

4 

taka  L  ■  9  and  N  ■  8  or  larger  to  ensure  retention  of  all  terms  0(a  )  . 

5.  Stability 

We  wish  to  consider  the  stability  to  infinitesimal  two  dimensional 
disturbances  of  the  steady  shapes  calculated  by  the  method  described  in  §  4. 

It  is  assumed  that  the  boundaries  are  deformed  by  disturbances  and  that  the 
vortidty  inside  the  vortices  remains  constant.  There  are  of  course  a 
continuum  of  oscillations  associated  with  changes  of  the  vortlclty 
distribution,  but  these  are  not  relevant  to  the  intrinsic  stability  of  the 
street  to  disturbances  produced  by  motion  of  boundaries  or  the  action  of 
conservative  forces.  We  suppose  now  that  the  coefficients  of  the  expansions 
(4.1),  etc.,  are  perturbed  by  Infinitesimal  functions  of  time  which  are 
denoted  by  primes.  The  linearised  unsteady  boundary  condition  that  the  vortex 
boundaries  are  material  boundaries  gives  the  equation 

dz1*  dz_  .  dz*  dz* 
dt  C  dt  dt  5  d? 

-  (u*  -  iv*)c  £  +  (u‘  +  iv*)  ?  ^ 

+  (u  -  iv)  c  +  (u  +  lv)  5  4z^  (5.d 
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Co  be  satisfied  on  each  vortex  of  Che  rows*  Unprimed  quantities  denote  the 
undisturbed  steady  state. 


We  now  have  to  handle  the  problem  that  for  subharmonic  disturbances  the 
oscillations  of  each  vortex  may  be  different,  so  that  we  appear  to  be  faced 
with  infinitely  many  equations  for  Infinitely  many  unknowns.  Irrespective  of 
any  truncation  that  is  employed.  This  difficulty  is  overcome  as  follows.  The 
velocity  Induced  by  a  vortex  is  an  analytic  function  of  the  coefficients. 

Let  (u  -  lv)  denote  the  velocity  Induced  by  the  vortex  centered  at  m & . 

n 

Then,  for  example,  the  change  in  this  velocity  at  a  fixed  value  of  z 
produced  by  the  perturbation  a'n>m  Is 


(u*  -  iv') 


d(u  ~  lv)  a, 

da  n,m 

n 


(5.2) 


where  u  -  iv  is  given  by  (4.5)  as  a  function  of  the  undisturbed 
coefficients.  Because  of  the  periodicity  of  the  undisturbed  motion,  the 
derivative  in  (5.2)  is  the  same  function  of  z  -  ml  for  every  vortex  in  the 
same  row.  Expressions  similar  to  (5.2)  hold  for  the  conjugate  velocity  and 
the  velocities  Induced  by  the  vortices  in  the  other  row.  It  follows  therefore 
that  normal  mode  solutions  of  (5.1)  can  be  obtained  by  restricting  attention 
to  one  vortex  in  each  row,  namely  the  vortices  given  by  (4.1)  and  (4.2),  and 
supposing  that  the  perturbations  of  the  other  vortices  are  related  to  the 
perturbations  of  these  vortices  by  relations 


a*  -  a*  elmf  ,a'*  -  a'*  e^.R'.  -  R'.e11"*,  R.’*  -  R’*  einf  (5.3) 

n,m  n  n,m  n  l,m  1  l,m  1 


b'n  »  "  bneiaf’  *  b'S  •** »  RS  a  "  ’  R2*m  “  R*2  ^  (5*4) 

n  j  id  n  D|0  n  }  m  z  z  |D  z 
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s 


where  f  is  related  to  the  subharmonic  wave  number  p  by 


f  -  2*p 


(5.5) 


Notice  the  subtle  point  that  although  the  undisturbed  solution  must  have 
the  property  that  aQ  and  a*Q,  etc.,  are  complex  conjugates  for  it  to  be 
physically  sensible,  this  need  not  be  true  of  Infinitesimal  perturbations. 
There  is  no  contradiction  here,  as  will  be  explained  in  more  detail  later.  I 
is  of  course  just  the  manifestation  in  the  present  context  of  the  fact  that 
the  eigenfunction  of  a  real  problem  may  be  complex. 

The  perturbed  velocity  u*  -  lv'  can  now  be  calculated.  We  have  the 
undisturbed  velocity  Induced  by  the  vortices  as  given  by  (4.12), 


-iai.  •  • 

.  -  lv  -  -J.  I  I 
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z’  -  R'1C(1  +  a./?  +  a-/c2+ 


1 


where 


(5.8) 


+  RjtCl  +  a' t/?  +  a’2/<;2+  ) 

and  g'n,  etc.,  are  the  perturbations  for  the  vortices  (4.1)  and  (4.2)  and 
are  obtained  from  (4.7),  etc.  The  differentiation  with  respect  to  z  in 
(5.7)  is  done  because  the  velocity  has  to  be  evaluated  on  |? j  ■  1  and  the 
boundary  condition  (5.1)  will  be  evaluated  in  terms  of 

The  boundary  condition  is  now  linear  in  all  the  perturbations  and 
solutions  are  searched  for  in  which  all  primed  quantities  are  proportional 
to  exp(ot).  If  we  denote  the  variables  R'i»  a’n,  R'*i,  a'*n,  R'2»  b*n  R'*2> 
b'*n  by  a  vector  the  boundary  condition  gives  an  eigenvalue  eq  latlon  of 
the  form 


a  H  (?)  -  M  (5.9) 

where,  at  least  in  principle,  the  functions  M  and  N  are  known  functions  of 
the  undisturbed  shape  and  the  subharmonic  wavenumber  p.  (The  infinite  sums 
with  respect  to  m  in  (5.7)  can  be  done  in  closed  form,  see  Appendix  B.) 

Again  in  principle,  the  eigenvalue  o  can  be  developed  as  a  series  in 
Inverse  powers  of  l ,  but  the  algebraic  complexity  is  fearsome  and  numerical 
methods  were  employed.  Corresponding  to  a  steady  solution  with  some  N  and 
L,  the  same  truncation  was  applied  to  the  perturbations,  i.e.  terms  with 
suffix  greater  than  N  were  discarded.  This  gives  a  column  vector  with 
4N  +  4  rows  for  The  forms  like  (4.15)  were  employed  for  the  calculation 
of  the  velocity  and  its  perturbation,  which  leads  to  an  expression  of  the 
boundary  condition  as  a  series  of  positive  and  negative  powers  of  ^  .  4N  +  2 
equations  are  obtained  by  equating  the  coefficients  of  ?n  to  zero  for  n  in 
the  range  -N  <  n  <  N  for  each  vortex.  Two  further  equations  are  obtained  by 
fixing  the  phase  by  requiring 


5 


The  problem  has  now  become  one  of  calculating  the  eigenvalues  of  a  generalised 
complex  eigenvalue  equation  of  order  4N  +  4,  which  can  be  approached  In 
standard  ways.  No  difficulty  was  found  in  accomplishing  this. 

Some  comments  are  in  order.  First,  It  will  be  noted  that  the  boundary 
conditions  for  the  unsteady  calculation  are  different  from  those  of  the  steady 
flow.  In  the  latter  case,  it  Is  necessary  to  ensure  that  the  equations  are 
Independent.  This  is  unnecessary  for  unsteady  calculations,  where  all  that  Is 
needed  Is  to  ensure  that  the  equations  are  not  Inconsistent.  Therefore  any 
Integral  constraint  that  leads  to  dependence  of  steady  equations  will  be 
automatically  satisfied  by  a  consistent  evolution  equation.  We  therefore 
retain  the  equations  that  come  from  the  terms  Independent  of  c  and  the 
coefficients  of  C  and  C for  the  second  vortex.  Satisfaction  of 
the  c  Independent  terms  In  equivalent  to  conserving  area  in  unsteady  flow,  and 
It  can  be  verified  from  the  computer  results  that  all  modes  with  a  *  0  keep 
the  area  constant.  Indeed,  this  provides  a  check  on  the  accuracy  of  the 
eigenvalues  and  eigenvectors.  As  described,  the  calculation 
keeps  and  <i>2  constant.  There  is  a  neighboring  equilibrium  state  In  which 
the  velocity  of  the  array  is  the  same  (Q  and  Q*  are  not  perturbed)  and  the 
strengths  and  sizes  of  the  vortices  are  altered.  We  expect  corresponding  to 
this  state  that  there  will  be  two  zero  eigenvalues  in  which  the  eigenfunction 
does  not  conserve  area.  These  are  indeed  found.  They  could  have  been 
eliminated  by  imposing  area  conservation  to  reduce  the  number  of  unknowns,  but 
this  would  have  been  less  convenient.  The  retention  of  (5.10)  gives  two 
Infinite  eigenvalues,  and  it  proved  easier  to  make  the  eigenvector  satisfy 
(5.10)  throughout  the  calculation  and  actually  solve  a  4N  +  2  order 


system*  This  gives  N  quartets  of  eigenvalues,  one  of  which  is  the 

cooperative  mode  and  the  others  are  the  first  N— 1  oscillation  modes,  and  the 

two  neutral  modes  in  which  the  area  changes . 

Secondly,  we  find  perturbations  which  contain  components 

like  a'  eot  +  and  eat  +  im*  for  the  m-th  vortex,  and  these  are  not 

n  n 

complex  conjugates  of  each  other.  To  avoid  confusion,  we  denote  the 
perturbation  to  a*n  by  •  They  produce  contributions  to  z'n  and  z'*m 
which  are  respectively 


(f)R.  .«<*>*  +  *"f 
n  l 


(5.11) 


3Wn’le°(f)t  +  iof 


(5.12) 


Now  change  f  into  -f  (or  2s  -  f).  We  get  the  contributions  to  z'm  and 
z'*m 

a'  (-f)R1eo(“f)t  ’  lof 
-  n--_ r  (5.13) 

C  1 

rnR*1?n-1ea(‘f)t  '  imf  (5.14) 


Because  of  the  symmetry  expressed  by  (2.2),  o(-f)  ■  o*(f),  and 

a'  (-f)  ■  (¥'  (f))*,  where  the  star  here  denotes  actual  complex  conjugate, 
n  n 

Thus  (5'.11)  is  the  complex  conjugate  of  (5.14),  and  (5.12)  is  the  complex 
conjugate  of  (5.13).  Physically  realistic  perturbations  are  therefore 
obtained  by  adding  (5.11)  and  (5.13)  for  the  disturbance  to  zm  and  (5.12) 
and  (5.14)  for  the  disturbance  to  zm  and  (5.12)  and  (5.14)  for  the 


disturbance  to  z*  . 


6 .  Results 


Steady  shapes  were  calculated  for  a  range  of  o  and  ic  using  the  tact  hod 

°*  $4  and  the  stability  of  these  shapes  to  inflnltesiaal  disturbances  was 
calculated  using  the  aethod  of  section  5.  For  the  majority  of  the 

calculations,  the  value  of  H  ■  8  and  L  •  9  were  chosen*  Results  were 
checked  in  two  ways.  First,  calculations  were  done  with  N  -  16,  L  -  17, 
and  N  ■  32,  L  -  33.  Second  an  alternative  code  which  employed  the 
'interior'  napping  of  the  exterior  of  the  vortex  into  the  interior  of  the  unit 
circle  was  employed  Instead  of  (4.1).  This  aethod  was  used  with 
L  -  N  -  12.  For  a  less  then  0.1,  there  was  no  significant  difference  between 
any  of  the  results.  Calculations  using  N  -  8  and  L  •  9  required  about  15 
minutes  of  CPU  for  the  computation  of  the  steady  shapes  and  about  the  same 
tine  for  the  computation  of  the  eigenvalues  on  a  VAX  11/750  computer.  Most  of 
the  tlae  went  for  the  calculation  of  the  Jacobian  and  the  matrices.  A  large 
amount  of  data  was  amassed.  Here  we  summarise  the  salient  features. 

The  behavior  of  the  shape  modes  showed  no  surprises.  In  fact,  the  effect 
of  finite  area  on  the  oscillation  frequencies  was  generally  insignificant. 

The  modes  were  all  stable  and  what  small  change  there  was  due  to  finite  area 
was  consistent  with  the  symmetry  requirements. 

The  Interest  lay  entirely  in  the  behavior  of  the  cooperative  modes  and  in 
particular,  the  effect  of  finite  area  on'  the  stability  boundaries  in 
the  p  -  k  plane  for  fixed  area.  A  typical  set  of  results  is  shown  in  figure 
7.  These  are  for  the  case  a  -  0.05.  The  ratio  r  ortex  radius  to 
longitudinal  separation  Is  0.1262.  The  solid  lines  show  the  stability 
boundary  in  the  vicinity  of  the  critical  point.  The  dashed  lines  are  the 
stability  boundary  according  to  equation  (3.30)  of  Rida  (1982).  It  will  be 
seen  that  the  numerical  results  indicate  that  the  stability  boundary  is  the 


Figure  7 . 

Numerical  results  for  the  stability  boundary  In  the  »c-p  plane  for  a  -  0.05.  Soli 
line  Is  result  of  calculation.  Dashed  line  is  Kida’s  prediction.  Stable  regions 
are  above  and  below  the  cross.  Also  shown  Is  the  cross  for 


zero  a  . 


degenerate  case  of  figures  5  and  6,  In  which  the  hyperbola  is  a  pair  of 
straight  lines  seating  at  the  center  px,  k^.  In  this  case,  the  dependence 
of  growth  rate  a  on  p  for  k  -  would  be  the  degenerate  case  of  figure  3  in 
which  the  ovals  shrink  to  two  points  at.  py  and  1-py.  The  existence  of  the 
cross  is  confined  by  the  plot  in  figure  8  of  the  difference  in  the  values 
of  p  on  the  stability  boundary  versus  k  very  close  to  py  and  k^.  The 

numerical  points  wars  found  to  lie  on  two  straight  lines  which  net  at  a  value 
of  the  difference  indistinguishable  froa  sero.  Similar  results  were  found  for 
all  other  values  of  a  employed.  In  fact,  some  runs  were  done  for  the 
relatively  large  value  of  a  ■  0.1854  (using  N  -  32),  and  showed  the  same 
behavior.  These  results  imply  that  the  street  is  always  unstable  except  for 
the  special  value  k  ■  Ky. 

For  the  values  of  k  shown  in  figure  8,  Kida's  results  imply  that  the 
street  la  stable.  (Note  that  the  results  drawn  in  figure  7  show  that  the 
pairing  mode  with  p  -  0.5  is  stabilised  by  finite  area,  as  calculated  by 
Saffnan  and  Schatsman  1982).  This  disagreement  is  not  unexpected,  since  as 

A 

argued  in  section  3  Kida's  calculation  is  not  consistent  to  0(a  )  as  is 
necessary  to  determine  the  behavior  near  the  center  of  the  hyperbolae. 

However,  according  to  the  argument  summarized  by  equation  (3.18),  the  change 
in  the  stability  boundary  away  from  the  cross  would  be  of  order  a  and  should 
be  correctly  calculated  to  this  order  by  Kida's  theory.  The  discrepancy 
between  the  numerical  results  and  Kida's  results  shown  in  figure  7  is  due  to 
the  value  of  a  being  too  large  for  the  perturbation  theory  to  be  accurate. 

This  is  confirmed  by  the  results  for  a  "  0.025  shown  in  figure  9,  where  the 
disagreement  is  confined  to  the  vicinity  of  the  cross. 

The  arguments  of  section  3  imply  that  px  -  1/2  and  fy  “  <c  should  be 
2 

proportional  to  o  •  Results  are  given  in  table  1  and  shown  in  figure  10.  It 
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Figure  9. 

Stability  boundary  for  a  ■  0.025.  Klda's  predictions  now  agree  with  the  nuoerical 
results  away  froa  the  cross  • 
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Figure  10. 


Foal don  of  Che  center  of  the  cross  for  varloua  a .  Solid  llnea  show 
_ i  1/2  .  ,1/2 

Pj  and  Kj  versus  a.  Dashed  line  Is  tangent  at  origin. 


vat  found  that  p'j  was  closer  to  quadratic  than  ic  ^  which  deviated 
considerably  for  the  larger  a.  Assuming  to  a  first  approximation 
that  ic'g  and  p'g  are  quadratic  plus  cubic  in  a,  we  have  from  the  results 
for  a  ■  0.005  and  a  •  0.0125  that  X  ■  0.557  and  y  ■  4.63.  It  follows  from 
equations  (3.14)  and  (3.15)  that 


C3  ■  13.6,  c$  ■  6.49, 


(6.1) 


The  agreement  of  these  values  with  those  calculated  by  Klda  (see  equation 
3.11)  is  very  good.  However,  the  existence  of  the  crosses  Implies  that 
D  -  0,  and  hence  from  equation  (3.13) 


-  -11.2 


(6.2) 


The  value  of  C5  agrees  well  with  the  value  predicted  from  the  Saffman  and 

Schatzman  results  for  p  -  1/2.  The  Important  result  is  the  difference 

between  (6.2)  and  Klda's  value  given  In  (3.11),  which  leads  him  to  predict 

stabilization  In  contrast  to  the  numerical  calculations.  A  consistent 

4 

perturbation  expansion  to  order  a  is  of  course  required,  but  this  Is  no 
simple  matter. 

It  may  be  of  interest  to  give  a  trivial  example  which  Illustrates  the 
consistency  problem.  Suppose  we  have  the  quadratic  equation 


+  ax  +  b  -  0, 


a  -  -2  +  E  +  0  (e2),  b  -  1  -  e  +  0(e2) 


Then  x  *  1  +  0(e),  but  use  of  the  quadratic  formula  shows  immediately  that 

2 

the  0(e)  term  cannot  be  calculated  without  knowledge  of  the  0(e  )  terms  In 
a  and  b. 

A  strange  property  of  our  results  is  that  there  was  no  sign  for  any  value 
of  a  of  the  crosses  breaking  up  into  hyperbolae.  This  implies  that  the 
functions  Hj  and  H3  are  for  finite  area  of  the  form 


H(p^c^x)  -  Dt(x  -  *x)2  -  D2(p  -  px)2  (6.3) 


where  Dj  and  O2  are  functions  of  p,K,a,  and  Kg,  px  are  functions 
of  a  only.  A  general  result  of  this  kind  should  be  capable  of  simple  proof, 
and  not  require  extensive  computation. 

It  is  known  experimentally  that  the  observed  values  of  k  in  the  vortex 
street  wake  of  a  cylinder  increase  with  downstream  distance,  as  presumably 
does  the  area  due  to  viscous  diffusion.  It  would  be  of  Interest  to  determine 


if  the  variation  of  k  and  a  were  such  as  to  keep  the  street  in  the  state  of 
stability  as  described  by  the  plot  in  figure  10,  but  this  is  at  present 
difficult  to  do  as  the  sice  of  the  vortices  is  not  easily  estimated  from  the 


data. 
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Appendix  A  Kinematic  Constraints 

Consider  the  double  row  with  vortices  of  strengths  and  rj  in  two 
parallel  horizontal  rows*  Suppose  that  the  velocity  at  Infinity  is  U^,V 
above  the  top  row  and  U2,V  below  the  bottom  row.  Conservation  of  mass 
requires  that  V  is  the  same  above  and  below.  The  equivalence  of  circulation 
and  vortlclty  flux  gives  the  relation 

u2  -  Di  -  (r1  +  r2)/i  (A1) 

where  l  is  the  distance  between  neighboring  vortices  In  the  same  row.  We 
note  the  vector  Identity 


u  x«  -  ZlrH2)  *  I»  (A2) 

Further,  for  two-dimensional  flow,  •  uik_,  where  jc^  Is  normal  to  the  plane. 
Then  using  Green's  theorem. 


/  u  x  u  dxdy  ■  -k  x  /u  u>  dxdy 

1  2 

■/  j  i  »  ds  -  /u  (u  *  n)ds 

-  1  V  l  <U2  -  Ui>  1  T  (Ui2  -  U22)  <A3> 

where  the  integral  is  over  one  wavelength  of  the  flow. 

Now  the  left  hand  side  of  (A3)  is  proportional  in  magnitude  to  the 


velocity  of  the  vortlclty  centroid,  and  this  is  zero  if  the  flow  is  steady. 
Hence  in  this  case,  either 


ox  -  -u2.  V  -  0, 


o.  -  o,  ,  r.  -  -r,. 


The  former  cat*  is  the  mixing  layer,  where  the  speed  of  the  vortices  is  the 
arithmetic  wean  of  the  speeds  on  the  two  sides.  The  latter  is  the  wake,  where 
the  vortex  speed  is  not  the  arithmetic  aean  and  the  circulations  are  equal  and 
opposite.  These  results  show  that  there  cannot  be  a  continuous  family  of 
steady  solutions  going  from  wake  to  mixing  layer  with  U2  changing 
contlnuoualy  from  U}  to  -U^ . 

Of  importance  for  our  numerical  method  is  the  fact  that  for  the  wake 
where  (A5)  holds,  the  right  hand  side  of  the  identity  (A3)  vanishes  and  hence 
generally,  writing  separately  the  contributions  to  the  left  hand  side  from  the 
vortices  in  each  row,  we  have 


J  u.  «  dxdy  +2 /  u  u  dxdy  ■  0 


which  is  a  consequence  of  the  kinematic  relation  between  velocity  and 
vortlclty.  Thus  whenever  the  vorticity  field  is  chosen  so  that  one  part  of 
(A6)  vanishes,  as  is  required  for  the  solution  to  be  steady,  the  other  side 
must  also  vanish  automatically.  Hence  equating  separately  to  zero  the  parts 
of  (A6)  in  a  calculation  of  steady  flow  does  not  give  independent  equations . 
This  is  the  Justification  for  dropping  the  equations  for  the  coefficients 
of  c  and  (  *  on  the  second  row  vortices. 


Appendix  B  Sun  Formulae 

The  calculation  of  the  steady  shapes  and  stability  requires  the 
evaluation  in  closed  forms  of  the  Infinite  suma 


l 

m*-«»  (tn+e  ) 


where  s  is  an  arbitrary  complex  number,  f  is  an  arbitrary  real  number 
and  n  is  a  positive  integer.  If  s  ■  0,  the  term  with  m  ■  0  in  the  sum 
is  excluded. 

It  is  convenient  to  deal  separately  with  s  ■  0  and  s  *  0.  For  the 
first  case,  we  have  for  n  ■  1,  by  direct  summation  or  the  elementary  theory 
of  Fourier  series, 


-  i(ir -f )  0  <  f  <  2ir 


■  0  for  f  -  0  or  f  -  2*  . 


Values  outside  this  range  are  obtained  by  applying  periodicity  in  f. 

For  n  >  l,  we  integrate  (B2)  n  -  1  times  with  respect  to  f,  choosing  the 
arbitrary  constants  so  that  the  resulting  expression  is  2ir  periodic  in  f. 
For  s  *  0,  we  have  the  formulae 


imf  is(2ir-f) 

(m+8 )  ^  *  2rr  is  T" 

e  -  1 


0  <  f  <  2ir  , 


■  w  cot  irs  f  ■  0  or  f  ■  2ir  . 


Differentiation  n  -  1  times  with  respect  to  s  gives  the  sums  for 
arbitrary  n.  Note  that  it  is  only  for  n  -  1  that  the  sum  is  not  a 
continuous  function  of  f • 

The  processes  of  integration  and  differentiation  are  easily  automated  and 
can  be  evaluated  by  symbolic  manipulation  programs • 
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20.  ABSTRACT  (Continue  on  rereree  tide  II  neeeeeery  and  Identity  by  block  number) 

The  stability  of  two-dimensional  infinitesimal  disturbances  of  the  inviscid 
Karman  vortex  street  of  finite  area  vortices  is  reexamined.  Numerical  results 
are  obtained  for  the  growth  rate  and  oscillation  frequencies  of  disturbances 
of  arbitrary  subharmonic  wavenumber  and  the  stability  boundaries  are 
calculated.  The  stabilization  of  the  pairing  instability  by  finite  area 
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20.  ABSTRACT  (cont.) 


demonstrated  by  Saffman  and  Schatzman  (1982)  is  confirmed  and  also  Kida’s 
(1982)  result  that  this  is  not  the  most  unstable  disturbance  when  the  area  is 
finite.  Contrary,  however,  to  Kida's  quantitative  predictions,  it  is  now 
found  that  finite  area  does  not  stabilize  the  street  to  infinitesimal  two- 
dimensional  disturbances  of  arbitrary  wavelength  and  that  it  is  always 
unstable  except  for  one  isolated  value  of  the  aspect  ratio  which  depends  upon 
the  size  of  the  vortices. 


